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SUMMARY 


This  report  presents  a  partial  description  and  a  discussion  of  the  quantitative  results  based 
on  the  exact  solution  of  the  following  theoretical  problem:  For  a  sphere  of  arbitrary  size, 
composed  of  homogeneous  material  with  a  finite  dielectric  constant  and  conductivity,  receiv¬ 
ing  a  constant  flux  of  energy  from  a  single  direction  in  the  form  of  plane  electromagnetic 
waves  of  a  given  frequency  and  state  of  polarization,  find  the  total  amount  of  energy  absorbed 
and  scattered  by  the  sphere,  as  well  as  the  specific  intensity  and  the  state  of  polarization  of 
the  energy  scattered  in  a  given  direction,  at  a  large  distance  from  the  sphere.  A  complete 
arulytical  solution  to  this  problem  was  obtained  some  time  ago  by  Gustav  Mie  on  the  basis 
of  Maxwellian  field  theory.  The  present  study  gives  accurate  numerical  results,  based  on 
Mie*s  expressions,  for  a  wide  range  of  basic  parameters  not  treated  earlier.  Illustrated  and 
discussed  by  means  of  selected  examples  are  the  effects  of  variations  in  dielectric  and  con¬ 
ducting  properties  on  the  total  scattering  and  absorption  cross  sections  as  a  function  of  the 
relative  size  of  the  sphere,  as  well  as  on  the  differential  amplitude,  intensity,  and  polariza¬ 
tion  of  the  scattered  energy. 

The  work  was  undertaken  as  a  preliminary  and  necessary  step  in  a  general  study  of 
planetary  atmospheres  containing  particles  of  various  sizes  and  types,  whose  characteristics 
are  to  be  determined  by  means  of  the  reflected  and  transmitted  visible  and  infrared  sun¬ 
light.  The  results,  after  proper  scaling,  are  equally  applicable  to  spheres  receiving  electro¬ 
magnetic  radiations  of  much  lower  frequencies,  such  as  microwaves. 
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The  main  results  of  this  work  have  appeared  in  D.  Deirmendjian,  R.  Clasen,  and  W.  Viezee, 
"Mie  Scattering  with  Complex  Index  of  Refraction,"  /.  Opt,  Soc.  Am,,  Vol.  31,  No.  6,  June, 
1961,  pp.  620-633.  The  minor  errors  in  the  curves  in  Figs.  6,  7,  9,  and  10,  noted  in  that  pub¬ 
lication,  have  been  corrected  in  the  corresponding  figures  of  the  present  report. 
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INTRODUCTION 


In  recent  years  an  increasing  number  of  computations  and  tables  related  to  the  scattering 
of  light  on  homogeneous  dielectric  spheres  have  appeared  in  the  literature.  A  comprehensive 
list  and  an  evaluation  of  these  appear  in  van  de  Hulst's  authoritative  monograph/'^  Since 
the  publication  of  his  study  in  1937,  additional  work  has  been  done,  but  in  most  of  the 
published  material,  only  the  total  scattering  cross  section  is  listed,  together  with  individual 
Mie  coefficients,  for  the  case  of  nonabsorbing  spheres  of  various  sizes  and  indices  of  refrac¬ 
tion.  The  amplitude  and  intensity  of  the  scattered  light  as  a  function  of  scattering  angle 
are  usually  omitted.^ 

A  complete  knowledge  of  the  intensity  and  polarization  of  light  that  has  been  scattered 
in  various  directions  is  essential  in  understanding  the  role  of  large  particles,  such  as  aerosol 
and  cloud  particles,  in  the  radiative  transfer  of  visible  and  infrared  light  in  atmospheres 
containing  such  particles.  Moreover,  a  consideration  of  dielectric  spheres  alone  is  not  suffi¬ 
cient.  For  example,  it  has  been  shown*’’ that  in  the  scattering  of  infrared  radiation  by  cloud 
particles,  a  complex  index  of  refraction  is  appropriate,  since  water  and  ice  display  infrared 
absorption  bands.  An  understanding  of  the  scattering  and  absorption  of  ultraviolet,  visible, 
and  infrared  light  by  partially  absorbing  particles  is  also  important  in  problems  presented 
by  sunlit  planetary  atmospheres  other  than  our  own,  as  well  as  by  interplanetary  and  inter¬ 
stellar  particles  of  a  size  comparable  to  the  wavelength.*^’ 

We  might  sketch  the  extent  of  present  knowledge  of  the  scattering  mechanism  (more 
fully  desaibed  by  van  de  Hulst,*^’  Twersky,*”  and  others)  as  follows;  The  particles,  assumed 
spherical,  are  characterized  by  two  parameters,  x  and  m.  The  size  parameter,  x,  is  the  product 
of  the  free-space  propagation  constant,  k,  and  the  radius,  4,  of  the  sphere.  The  complex  index 
of  refraction  m  is  related  to  the  dielectric  constant  and  to  the  conductivity  of  the  substance. 
Using  these  two  parameters,  we  may  classify  the  types  of  scattering  into  the  following  prin¬ 
cipal  cases. 

tTbese  ctmarkt  do  not  apply  to  the  Ittertture  appearing  after  March,  1961. 
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UGHT  SCATTERING  ON  SPHERES  OP  FINITE  SIZE 


RAYLEIGH  SCATTERING 

Rayleigh  scattering  is  restricted  to  the  case  x  1,  but  m  is  arbitrary.  This  limiting 
case,  applied  to  a  purely  molecular  atmosphere,  successfully  explains  most  observed  features 
of  the  color,  the  intensity,  and  the  polarization  of  the  sunlit  sky.  It  is  also  applicable  to  prob¬ 
lems  of  radar  backscatter  from  cloud  and  rain  droplets,  provided  their  dimensions  are  suffi¬ 
ciently  small  compared  with  the  wavelength. 


RAYLEIGH-GANS  SCATTERING 

In  Rayleigh-Gans  scattering  we  have  I  w  —  1  I  1  and  2x  |  m  —  1 1  1.  This  is 

a  limited  extension  of  the  Rayleigh  case.  The  conditions  are  more  stringent,  in  that  the  prop¬ 
erties  of  the  sphere  must  differ  only  slightly  from  those  of  the  surrounding  medium.  The 
restrictions  on  size,  however,  are  less  stringent  than  in  the  Rayleigh  scattering  described  in 
the  preceding  section,  since  the  case  x  >  1  can  be  treated,  provided  the  second  condition 
above  is  satisfied.  The  approach  is  useful  in  treating  the  radar  backscatter  from  cloud  and 
rain.  Furthermore,  such  nonspherical  particles  as  cylinders  and  elipsoids  can  be  considered, 
where  the  radius  of  a  sphere  is  replaced  by  another  characteristic  dimension. 


VAN  DE  HULST  APPROXIMATION 

The  van  de  Hulst  approximation  is  restricted  to  |  w  —  1 1  <<  1,  but  x  may  be  arbitrary. 
This  approximation  has  proved  most  successful  in  predicting  the  total  extinction  cross  section 
and  its  dependence  on  the  size  of  the  spheres.  The  theory  is  less  successful  in  the  case  of 
absorbing  spheres.  Essentially,  the  same  type  of  approximation  is  involved  in  what  is  known 
as  the  W.K.B.  approach,  which  describes  the  angular  scattering  and  polarization  as  well, 
within  small  angles  of  scattering.  The  theory  and  applications  have  been  discussed  in  Refs. 
3  and  6  through  9.  Methods  of  improving  the  approximation  have  been  discussed  by  Sekera 
in  Ref.  10.  The  method  so  far  fails  to  give  a  correct  estimate  of  the  intensity  and  polarization 
at  large  angles  of  scattering,  particularly  in  the  backward  hemisphere. 

GEOMETRIC-OPTICS  LIMIT 

The  geometric-optics  case  corresponds  to  x  — >  oo  with  m  real  and  either  finite  or  in¬ 
finite.  This  case  is  equivalent  to  a  wavelength  that  is  infinitesimal  compared  with  the  size  of 
the  scatterer,  and  independent  rays  of  light  can  be  analy-zed  as  they  enter  and  leave  the 
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medium.  The  pattern  of  the  scattered  light  is  assumed  to  be  given  by  effects  of  refraction, 
reflection,  and  diffraction.  This  approach  has  provided  us  with  qualitative  explanations  of 
ordinary  rainbows,  diffraction  patterns,  and  halos,  produced  by  cloud  and  rain  droplets.  It 
fails  to  explain  the  total  cross  sections,  the  polarization  observed  in  the  rainbow,  and  other 
quantitative  details  in  the  scattered  pattern,  because  important  mechanisms,  such  as  contribu¬ 
tions  from  surface  waves  and  internal  losses  cannot  be  easily  considered. 

ELECTROMAGNETIC  THEORY — MIE  SCATTERING 

In  this  instance  both  m  (complex)  and  x  may  be  arbitrary.  For  the  case  of  spheres,  Mie 
solved  a  boundary-value  problem  using  Maxwell’s  electromagnetic  theory.  The  solution  is 
exact,  and  it  describes  the  field  both  inside  and  outside  a  sphere  for  any  scattering  angle, 
size,  and  index  of  refraction.  At  the  limits  x<^<[  1  and  x— >  oo,  it  gives  the  known  fea¬ 
tures  of  the  Rayleigh  and  geometric-optics  approximations.  Because  the  exact  theory  is  most 
useful  between  these  limits,  problems  in  this  range  are  usually  referred  to  as  Mie  scattering 
problems  and  the  corresponding  particles  as  Mie  particles.  This  is  the  theory  on  wliich  the 
present  computations  are  based. 

Calculations  of  Mie  scattering  with  complex  index,  especially  for  spheres  of  finite  rela¬ 
tive  size  X,  are  rather  scarce.  The  work  of  Schoenberg  and  Jung,  that  of  Schalen,  and  Mie's 
original  computations,  as  described  by  van  de  Hulst,^^^  generally  involve  small  spheres  (less 
than  one  wavelength  in  diameter),  so  that  a  series  expansion  in  x  may  be  used.*^^’  For  larger 
spheres,  we  have  a  few  limited  cases  given  in  Lowan’s  tables.*^*^  In  addition,  values  of  the 
extinction  cross  section  alone,  based  on  the  Mie  theory,  are  given  by  Johnson  and  Terrell,^*’ 
and  by  Haugen.^'®^ 

The  purpose  of  this  report  is  to  describe  new  and  more  detailed  computations,  per¬ 
formed  at  The  rand  Corporation,  of  the  scattering  and  absorption  characteristics  of  spheres 
of  various  optical  properties,  mainly  corresponding  to  liquid  water  illuminated  by  visible  and 
infrared  radiation.  Examples  are  presented  mostly  by  means  of  suitable  plots  of  the  functions. 


NOTATION  AND  COMPUTATIONAL  SCHEME 


The  mathematical  formalism  relating  to  the  Mie  theory  will  be  kept  to  a  minimum,  by 
referring  the  reader  to  van  de  Hulst's  clear  and  consistent  discussion.^  With  minor  varia< 
tions,  we  adopt  the  symbols  and  definitions  introduced  by  that  author.  In  brief,  we  are  given 
the  following:  the  size  parameter,  x  =  ka,  where  k  —  lir/X  and  X  =  the  free-space  wave¬ 
length  of  the  incident  radiation;  m  =  r  —  i/,  the  complex  index  of  refraction  in  the  interior 
of  the  sphere  with  respect  to  the  medium,  where  (e  —  47ria/it>y^^  e  being  the  dielectric 
constant,  a  the  conductivity,  and  w  the  circular  frequency  of  the  electromagnetic  wave;  and 
the  scattering  angle  B  measured  from  the  forward  direction. 

The  incident  plane-wave  radiation  is  assumed  unpolarized  so  that  the  electric  field  may 
be  represented  by  two  independent  and  mutually  perpendicular  oscillations  of  unit  ampli¬ 
tude  in  the  form 


propagating  in  the  direction  z.  The  resultant  field  will  be  the  sum  of  the  initial  field  plus 
a  scattered  field,  whose  electric  amplitude  may  be  expressed  in  the  form  of  two  vector  com¬ 
ponents,  A^  perpendicular  to  the  scattering  plane  (in  which  the  scattering  angle  6  is  meas¬ 
ured)  and  Aj  parallel  to  the  scattering  plane.  At  large  distances  from  the  sphere,  in  the 
so<allcd  radiation  limit,  these  vectors  have  no  component  along  the  direction  of  propaga¬ 
tion.  According  to  the  Mie  theory,***  the  magnitude  of  these  far-field  amplitude  vectors  may 
be  expressed  in  dimensionless  form  as  follows: 

kA,  =  5.(x.  w,  ^  {-.’r.  -1-  Krn) .  (1) 

€0  2  J 

~  »(» + 1)  *  (2) 

^Sce  pp.  of  Ref.  1. 
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and 


=  ^  £  (2«  +  1)  Re  {a„  +  b„), 

«■=! 

K.^{X,  (2«  +  0  {  I  ‘*n  r  +  I  1*  )■ 


where  5^  and  5,  are  the  dimensionless  complex  amplitude  fund 
total  extinction  and  scattering  cross  sections  divided  by  the  geotra^ 
the  ’’efficiency  factors”  and  2,0  in  van  de  Hulst’s  treatmecal 
are  the  Mie  complex  scattering  coefficients;  and  7r^(^)  and 
coefficients 

'^nip)=^Pnip). 


%(m)  =M’r,(M)  —  (1  —  p‘) 


where  ft  =  cos  ^  and  P»(m)  is  the  Legendre  polynomial  of  order 
The  most  troublesome  quantities  to  treat  computationally 
since  they  involve  spherical  Bessel  functions  of  real  and  coi 
the  corresponding  expressions  may  be  written  in  terms  of  ordinary^ 
kind,  of  orders  ;i  +  of  real  argument,  and  circular  and  hyperl 
using  the  theory  of  Bessel  functions, and  after  some  manipula. 
as  given  by  van  de  Hulst,^  we  have 


(3) 


(4) 


K„^  and  K„  are  the 
etrical  cross  section  (called 
<*,(»»,  at)  and  b^{m,x) 
are  the  angle-dependent 


1. 


Lxe  the  coefficients  a,  and 
iplex  arguments.  However, 
Bessel  functions  of  the  first 
•lie  functions.  For  example, 
ion  of  the  Mie  expressions 


+  j]  [/..h(^)  +  ]  -  Jr^Siixy 

/■(— )7-,»V4  W 

where  involves  the  Bessel  fimctions  with  complex  argument  —  and  may  be  reduced  to 

the  form 


A.(mx)  - - h 

'  mx  ' 


1n.%(,rnxY 


The  expression  for  b^  is  exactly  as  in  (5)  except  that  A^/tn 
it  occurs. 


As  replaced  by  mA^  wherever 


tSeep.  12}  of  Ref.  1. 
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Actually  the  expression  for  A^{mx)  takes  on  the  following  very  simple  recursion  form: 


A„{mx) 


— •  H - 

mx  n 
rnx 


(6) 


and  it  can  be  evaluated  for  any  positive  integer  n  from  the  tormula 


A„(mx) 


=:  ctn  (wx) . 


From  formula  (6)  it  is  clear  that  is  real  when  m  is  real,  and  (5)  reduces  to  the  cor¬ 
responding  expression  for  dielectric  spheres. 

The  half-order  Bessel  functions  of  argument  x  can  also  be  generated  with  the  usual 
recursion  formulas, and  the  problem  of  determining  the  real  and  imaginary  parts  of  any 
A^  is  one  of  rationalization  and  algebra.  The  coefficient  may  be  treated  similarly.  Further, 
the  angular  functions  and  reduce  to  Legendre  polynonuals  and  their  first  and  second 
derivatives,  which  again  are  amenable  to  recursion  techniques  and  well  adapted  to  digital 
computers.  (See  Appendix  A  for  details.) 

The  above  calculations  were  programmed  for  the  IBM-704  calculator  at  rand  in  the 
FORTRAN  language.  The  program  has  now  been  adapted  for  the  IBM-7090.  All  the  neces¬ 
sary  functions  were  generated  by  the  corresponding  recursion  formulas.  The  quantities 
and  were  terminated  cither  when 


Lf*]  10*^*,  or  when  n  =  1.2x  H-  9* 

n 

The  latter  criterion  is  somewhat  more  stringent  than  that  used  by  Penndorf  and  Goldberg/"^ 
who  have  published  detailed  tables  pertaining  to  the  coefficients  for  dielectric  spheres.  The 
values  from  the  IBM-704  checked  exactly  with  the  tables  cited  above  in  the  case  where  an 
identical  index  of  refraction  was  used.  Thus,  an  accuracy  of  at  least  five  decimal  places  was 
achieved  in  the  functions  4^,  x^/,  and 

In  the  choice  of  refractive  index,  we  were  guided  by  indices  summarized  by  Centeno^**^ 
for  liquid  water  in  the  visible  and  infrared  regions,  and  by  van  de  Hulst*'^  for  ir<Mi  sfdieres 
in  the  visible.  The  actual  indices  used  initially  are  listed  in  Appendix  B.  The  size  interval 
and  size  range,  as  well  as  the  interval  in  the  scattering  angle  0,  were  varied  to  obtain  a  good 
description  without  either  too  much  or  too  little  detail  in  the  numerical  results.  The  choice 
here  is  also  governed  by  the  numerical  integration  with  respect  to  particle  size  eventually 
required  to  obtain  the  resultant  parameters  for  a  sample  distribution  in  a  polydispersed  medium. 
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As  for  computing  time,  one  of  the  longest  runs  on  the  IBM-704,  involving  x  in  50(0.5)70 
and  6  in  0°(0.5)5(1)10(5)180°,  took  less  than  thirty  minutes.  (Note  that  for  x  =  70  about 
100  terms  in  the  series  appearing  in  expressions  (1)  to  (4)  have  to  be  evaluated  before  con¬ 
vergence  is  achieved.) 

Once  the  primary  quantities  (1),  (2),  (3),  (4)  are  known,  the  following  derived  quan¬ 
tities  may  easily  be  determined: 


ri(x,  w, «)  =  I S,  1’  =  SX .  (7) 

/,(x.  w,e)=  |5j*  =  5,5t,  (8) 

/3(x,w,fl)  =  Re{S,St},  (9) 

/,(x,  m,  ^)  =  ~Im  {SXi}^  (10) 

where  the  asterisk  denotes  the  complex  conjugate  of  the  quantity;  and 

^.b.(^.  w)  ““  (11) 


Here  and  are  proportional  to  the  intensities  of  the  light  scattered  per  unit  solid  angle 
into  the  direction  0  with  respect  to  the  direction  of  the  incident  light,  which  consists  of  two 
units  of  parallel  flux  of  neutral  radiation;  these  quantities  give  the  degree  of  linear  polariza¬ 
tion;  /g  and  are  related  to  the  Stokes  parameters  that  describe  the  plane  of  polarization  and 
the  ellipticity;  and  is  the  relative  absorption  cross  section  in  the  case  of  absorbing  particles. 

All  the  functions  (1)  to  (4)  and  (7)  to  (11)  were  printed  out  on  standard  IBM  sheets 
for  the  various  indices  of  refraction  and  ranges  of  parameters  shown  in  Appendix  B.  For 
obvious  reasons,  it  is  not  possible  to  present  here  a  complete  table  of  values.  In  the  follow¬ 
ing  sections,  we  shall  attempt  to  describe  various  features  of  the  results  by  means  of  selected 
examples.  For  the  substance  and  wavelength  corresponding  to  these  examples,  the  reader  is 
referred  to  Appendix  B. 
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THE  COEFFICIENTS  AND  K 


As  VAN  DE  Hulst  hu  showti/  the  coefficients  and  h,,  when  plotted  in  the  complex  plane, 
always  fall  on  a  circle  of  radius  0.)  with  the  center  at  (0.3, 0.0)  provided  m  is  real.  He  also 
pointed  out  that  this  is  no  longer  true  when  m  is  complex  and  the  absorption  is  large.  The 
present  computations  give  a  more  detailed  description  of  the  Mie  coefficients  for  partially 
absorbing  spheres. 

Figures  la  and  ih  show  the  continuous  variation  of  a^{x)  and  ^,(x)  as  a  function  of  x 
in  the  complex  plane  for  luxiabsorbing  and  for  slightly  absorbing  spheres.  In  the  notuibsorb- 
ing  case,  both  coefficients  trace  the  circle  clockwise.  When  a  small  amount  of  absorption  is 
introduced,  the  curves  begin  to  spiral  inward  from  the  circle,  and  a  series  of  small,  counter¬ 
clockwise  loops  appear,  which  increase  in  size  and  show  a  tendency  to  converge  toward  an 
area  near  the  center  of  the  circle. 

An  increase  in  absorption  index  from  0.0472  to  0.0643  (Fig.  1^)  has  the  effect  of 
amplifying  the  loops  and  "accelerating"  the  convergence  toward  an  area  near  the  center  of 
the  circle.  In  both  cases  ^i(x)  develops  its  loops  for  smaller  x  than  tf,(x);  otherwise  its  be¬ 
havior  in  the  complex  plane  is  quite  similar  to  4,(x).  The  higher  order  coefficients  also 
show  a  similar  behavior  when  plotted  in  this  way. 

The  case  of  very  large  absorption  (metallic  spheres)  has  been  adequatdy  described  by 
van  de  Hulst.*  Our  comparable  results  (not  shown  in  a  separate  diagram)  corroborate  die 
fact  that  tf,(x)  traces  a  single  loop,  and  thereafter  its  locus  rapidly  approaches  a  counter¬ 
clockwise  circle  u  x  increases.  In  the  case  of  /',(x),  the  same  limiting  circle  appears,  cen¬ 
tered  at  (0.3,  0.0)  and  with  a  radius  equal  to  one  half  the  absolute  value  of  the  Fresnd 
reflection  coefficient  for  perpendicular  incideiKe. 


m-l  ^r(f- 1)’ +  /’■]'* 

m  +  l  L(»-+l)‘  +  r'*J  ■ 


(12) 


tSeepp.  13).  276/.  6(M.  I. 

tSce  p.  276,  Fi/s.  33,  •«!  p.  279  of  Ref.  1. 
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Fig.  1 — ^The  effect  of  absorption  on  the  complex  coefficients  4i(x)  and  the  first  terms  in 

the  Mie  series.  G>mputed  points  are  represented  by  dots  and  corresponding  x-values  are  shown 
along  the  curves,  (a)  m  =  1.29  —  0.0472/  in  the  range  x  =  0(0.23)12.0.  The  outer  circle  corre¬ 
sponds  to  4i(x)  for  m  =  1.29  (no  absorption),  (b)  m  =  1.29  —  0.0643/  in  the  range 
X  ^  0(0.3)18.0.  The  outer  circle  corresponds  to  ^i(x)  for  m  =  1.29. 


Fig.  l — contmH0d 


From  the  foregoing,  it  may  be  presumed  that,  for  small  and  moderate  absorption  and  for 
sufficiently  large  values  of  x,  the  loops  will  also  degenerate  into  a  circle  with  diameter  given 
by  (12). 

Note  in  Fig.  1  that  for  a  given  x  the  points  <»,(x)  and  b^{x)  for  m  =  r,  and  m  =  r  —  // 
fall  on  the  same  radius  drawn  from  the  center  of  the  circle.  This  indicates  that  the  absotp- 
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tion  reduces  the  magnitude  of  both  I  a„  —  %  1  and  \  h^  —  |  by  a  factor  that  increases  with 

//r  and  with  x,  provided  y* / f  is  small.  The  damping  effect  of  the  absorption  index  /  on 
and  does  not  apply  for  all  x,  however.  For  sufficiently  small  x,  the  addition  of  absorp¬ 
tion  results  in  higher  values  of  these  coefficients,  particularly  b^. 

Another  way  of  plotting  these  coefficients  is  to  keep  x  constant  and  vary  ».  This  kind 
of  representation  (not  shown  in  a  separate  diagram)  again  shows  the  damping  effect  de¬ 
scribed  above,  the  effect  increasing  with  r//  but  not  necessarily  with  «.  As  n  increases,  the 
values  of  and  always  tend  to  zero  for  some  n  that  depends  on  x,  indicating  the  point  of 
convergence  of  the  series  in  expressions  (1)  to  (4).  In  general,  mterpolations  of  individual 
values  of  these  coefficients,  either  with  respect  to  x  and  or  r  and  /,  should  not  be  neces¬ 
sary  once  a  correct  computing  program  has  been  worked  out.  Interpolations  of  the  functions 
involving  the  summed  series  may  be  necessary  and  are  justified  by  the  smoothness  of  these 
functions. 


THE  AMPLITUDE  AND  INTENSITY  IN  THE  FORWARD  DIRECTION 

Figure  2  represents  a  smooth  curve  in  the  complex  plane  fitted  to  the  computed  values 
of  x"*5(0°)  as  a  function  of  x,  where  S(0°)  is  the  amplitude  of  the  wave  scattered  in  the 
forward  direction  (^  —  0°)  for  absorbing  and  nonabsorbing  spheres,  respectively.  In  this 
direction,  the  polarization  of  the  incident  light  is  preserved  for  all  complex  indices  of  refrac¬ 
tion.  The  case  m  =  1.29  may  be  compared  directly  with  that  for  m  =  1.33  shown  by  van  de 
Hulst.^  In  Fig.  2  the  values  extend  to  x=  22.3  with  intervals  of  0.23  in  x.  Loops  appear  for 
X  >  8  increasing  in  frequency  with  x.  For  x  — >  oo,  these  should  degenerate  into  the  point 

0),  which  is  the  geometric-optics  (or  high-frequency)  limit  given  by  the  diffraction 
term.'^‘* 

The  curve  for  1.29  —  0.0643/  shows  that  the  addition  of  absorption  in  general 
damps  the  amplitude,  that  the  loops  appear  at  larger  values  of  x,  and  that  the  diffraaion 
limit  is  reached  for  smaller  values  of  x  than  in  the  case  of  pure  dielectric  spheres.  The 
smoothness  of  this  spiral  indicates  that  it  is  possible  to  extend  the  analytical  approximations 
of  van  de  Hulst*  for  !  w  —  1 1  0  to  include  absorbing  spheres.  When  I  w  —  1 1  is  finite, 

the  van  de  Hulst  expressions  may  be  considerably  improved  by  means  of  a  correction  term/*’ 
However,  this  possibilit)*  is  limited  only  to  the  forward  amplitude  and  the  total  extinction 
cross  section  (see  the  following  section). 

^See  p.  264  of  Ref.  1. 
tSec  p.  176  of  Ref.  1. 
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The  effect  of  increasing  the  real  part  r  of  the  index  of  refraction,  but  keeping  /  con¬ 
stant,  is  to  introduce  waves  and  loops  in  the  smooth  spiral.  The  geometric-optics  limit  is 
approached  more  slowly  than  in  the  case  of  r  near  unity.  However,  this  limit  has  always  the 
value  V2  regardless  of  the  optical  properties  of  the  sphere. 

Figure  3  is  a  plot  of  the  quantity  equivalent  to  the  square  of  the  absolute  value 

of  the  amplitude  shown  in  Fig.  2.  This  quantity  is  proportional  to  the  theoretical  intensity 
of  the  light  scattered  exactly  in  the  forward  direction,  which,  of  course,  cannot  be  distin¬ 
guished  from  the  transmitted  light  experimentally.  The  damping  effect  of  the  absorption  on 
both  the  large-  and  small-scale  oscillations  is  clearly  seen  by  comparing  the  curves  for 
m  —  1.29  (full  line)  and  m  =  1.29  —  0.0645/  (dashed  line).  Increasing  the  real  part  of  m 
(dash-dot  curve)  is  equivalent  to  shifting  the  first  resonance  maximum  toward  smaller  sizes. 


Fig.  3 — The  logsrithm  of  the  intensity  function  at  aero  Ksttering  angle,  /r*/(0^)  as  a  function  of  x  for 
the  three  indice*  i  refraction  shown.  Points  computed  at  intervab  of  0.25  in  x  have  been  joined  by 
straight  1....  .  ine  minor  OKillations  in  the  full  line  curve  correspond  to  the  OKillations  and  loops  of 
the  amplitude  spiral  in  Fig.  2. 
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For  sufficiently  large  x,  all  curves  in  this  type  of  representation,  regardless  of  m,  approach 
the  straight  line  .v"‘/(0°)  =  0.25,  which  is  the  geometric-optics  limit. 

The  curves  in  Fig.  3  are  also  related  to  the  infrared  aureole  near  the  sun’s  pro¬ 

duced  by  atmospheric  water  droplets  (see  Appendix  B).  The  scatter eu  intensity  for  a  given 
droplet  size  is  determined  on  the  curves  by  the  equivalent  x  value.  Furthermore,  the  area 
under  the  curves  in  Fig.  3  is  a  measure  of  the  aureole  intensity  at  the  sun*s  limb,  assuming 
a  particle  concentration  that  varies  as  a~*.  in  the  latter  case,  two  conclusions  are  immediately 
evident:  (1)  The  aureole  intensity  at  small  angles  must  increase  with  an  increase  in  the  upper 
limit  of  particle  size.  (2)  If  the  range  in  size  would  be  wide  enough,  it  becomes  difficult 
to  distinguish  absorbing  and  nonabsorbing  particles  by  measuring  the  aureole  intensity  only 
at  small  angles.  However,  if  the  particles  are  relatively  small  (x  <1 10),  and  still  show  the 
a~*  distribution,  the  aureole  intensity  next  to  the  limb  should  be  larger  for  nonabsorbing 
(w=:1.29)  than  for  absorbing  (w  =  1.29  —  0.0645/)  particles.  On  the  other  hand,  if  the 
real  part  of  the  index  is  increased  (m  =  1.525  —  0.0682/),  the  smaller  particles  become  more 
efficient  scatterers  in  the  forward  direction. 

THE  EXTINCTION  AND  ABSORPTION  CROSS  SECTIONS 

The  extinction  cross  section  is  connected  to  the  amplitude  of  the  wave  scattered  in  the 
forward  direction  according  to  the  cross-section  theorem,  which  in  the  present  notation  takes 
the  form 

^f..t(»».x)  =  ^Re{S(0)}.  (13) 

Figures  4a  to  4c  show  examples  of  the  extinction  and  absorption  cross  sections  plotted  as  a 
function  of  van  de  HulsFs  normalized  parameter  2  =  2x(f  —  1).  With  an  appropriate  re¬ 
labeling  of  the  z-axis,  the  graphs  may  b^  made  to  represent  variations  of  the  cross  sections 
with  wavelength  or  with  the  size  of  the  particles.  The  three  cases  shown  correspond  to  refrac¬ 
tive  indices  with  identical  real  parts  but  increasing  imaginary  parts  representing  liquid  water 
droplets  illuminated  by  infrared  radiations  of  different  wavelengths  (see  Appendix  B). 

The  full  lines  with  dots  in  Figs.  4a  to  4c  show  the  exact  Mie  values  computed  according 
to  (3),  (4),  and  (11).  It  is  known  from  previous  analyses^ that  in  this  type  of  repre¬ 
sentation  the  extinction  shows  a  principal  maximum  around  z  =  4  with  a  damped  periodic 
shape  for  larger  values  of  z.  With  the  inclusion  of  increased  amounts  of  absorption,  the 
principal  maximum  is  damped  until  it  disappears  (Fig.  4c)  and  reduces  to  a  plateau.  The 
larger  the  absorption,  the  smaller  the  value  of  z  at  which  the  asymptotic  value  2  as 
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Fig.  4 — The  logarithm  of  the  extinction  and  absorption  cross  sections  as  a  function 

of  the  normalized  size  parameter  2x{r  —  1).  Full  lines  with  dots  indicate  exact  Mie  values.  Dashed 
lines  correspond  to  the  equivalent  analytical  function  in  van  de  Hulst's  approximation  for  |  m  |  1, 

computed  for  the  same  index  of  refraction.  The  three  sets  of  curves  refer  to  a  refractive  index 
r  =  1.315  and  an  absorption  index  of  r'  =  0.0143  {a),  0.1370  (^),  and  0.4298  (c),  respectively. 


z  cx)  is  reached.  Further,  as  pointed  out  in  the  previous  section  (see  also  Fig.  2),  the 
secondary  oscillations  in  characteristic  of  nonabsorbing  spheres,  all  but  disappear  when 
absorption  is  introduced.  In  this  respect,  it  should  be  noted  that  the  extinction  curves  com¬ 
puted  by  Johnson  and  Terrell^'^  contain  some  errors,  which  must  be  attributed  to  an  insuffi¬ 
cient  accuracy  in  their  generating  functions.  For  example,  the  apparent  discontinuity  shown 
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(6) 

Fig.  4 — continued 


in  their  curve  for  »j  =  1.29(1  ““0.05/’)  between  x  =  6.6,  and  x  =  7  was  not  reproduced 
when  we  repeated  this  case  in  our  program.  Compare  the  following  values  of  in  Table  1. 
The  comparison  with  our  values,  which  carry  an  accuracy  of  five  decimal  places,  shows  that 
one  must  generate  the  basic  Mie  functions  to  several  significant  figures  to  obtain  sufficient 
accuracy  in  the  summations.  Except  for  such  details,  the  values  of  Johnson  and  Terrell^*’  are 
essentially  correct. 


Table  I 

X„t  ior  OT=  1.29(1  —0.05/) 


X 


Source 

1 

2 

4 

6 

6.6 

7 

17 

Johnson  and 
Tertell<«> 

0.261 

0.864 

2.28 

3.07 

J.19 

3.06 

2.40 

The  present  report 

0.2496 

0.8563 

2.275 

3.007 

3.062 

3.055 

2.319 

RESULTS 
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Fig,  4 — continued 


The  dashed  curves  in  Figs.  4a  to  4c  show  the  cross  sections  computed  for  the  same 
indices  by  means  of  van  de  Hulst's^  analytical  approximations  without  correction.^®^  The  gen- 
eral  shape  of  the  exact  curves  is  reproduced,  but  there  is  a  systematic  deviation  especially  in 
the  curves  for  These  deviations  in  the  extinction  curves  seem  to  vanish  as  z  — >  oo  but 
not  necessarily  in  the  absorption  curves.  A  useful  extension  of  these  approximations  has  been 
discussed  elsewhere/^^  and  the  correction  factors  proposed  are  found  to  check  against  the 
present  exact  computations  of  the  extinction  cross  section. 

VARIATION  OF  THE  COMPLEX  AMPLITUDE  WITH  SCATTERING  ANGLE 

The  angular  dependence  of  the  complex  amplitudes  5t(0)  and  S,(0)  for  any  given  x 
and  m  may  be  represented  by  continuous  curves  on  the  complex  plane.  The  usefulness  of  this 


f 


*■ 


tSee  pp.  179  and  181  of  Ref.  1. 
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representation  in  interpolations  of  the  amplitude  with  respect  to  0  is  limited  to  relatively 
small  sizes,  say  x  8.  For  our  purposes,  it  is  more  interesting  to  compare  the  amplitude 
curves  for  dielectric  and  for  partially  absorbing  spheres  of  an  identical  size  and  refractivity, 
since  this  has  not  been  done  previously. 

Figure  5  shows  the  amplitudes  for  particles  with  x  =  1  for  the  three  indices  shown. 
For  real  index,  the  amplitude  plots  for  5i(6)  and  S^{8)  are  practically  straight  lines 
almost  parallel  to  the  imaginary  axis.  There  is  very  little  variation  in  the  magnitude 


Fig.  5 — The  effect  of  increasing  absorption  on  the  amplitude  functions  5i(^)  (dashed  line)  and 
5,(^)  (full  line)  plotted  on  the  complex  plane,  in  the  case  of  a  small  sphere  (x  =  1).  Dots  indicate 
computed  points  at  10^  and  20^  intervals  in  $,  Those  for  0^,  90^,  and  180^  scattering  are  labeled. 
The  three  curves  correspond  to  f  =  1.31)  and  0,  0.1370,  and  0.4298,  respectively,  as  indicated. 
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of  5i(0)  with  9,  but  that  of  S^(9)  is  tninimhed  near  9  =  90^.  This  feature  is  similar 
to  that  obtained  in  the  Rayleigh  limit,  (x  <<  l),  where,  in  the  first  approximation,  we  have 
S^(0)  +  2)  andS,(0)  ==  cos  +  2),  respectively. 

The  other  curves  on  Fig.  5  show  that  the  amplitude  at  all  angles  increases  with  the  absorp¬ 
tion  index,  mainly  because  of  a  larger  real  part.  Further,  all  three  cases  arc  represented  by 
almost  straight  lines,  with  the  trivial  conclusion  that  interpolation  with  9  becomes  very  easy 
when  X  is  small.  For  any  x  1,  the  expansion  of  the  amplitude  in  power  series  of  x,  men¬ 
tioned  by  van  de  Hulst,^  and  used  by  Penndorf^^^  in  calculating  the  cross  sections,  may  be 
adequate  for  the  amplitudes  as  well. 

For  greater  values  of  the  size  parameter  x,  the  amplitudes  show  large  fluctuations  with 
9y  particularly  in  the  backward  hemisphere.  This  creates  a  problem  in  the  choice  of  an  interval 
L9  small  enough  to  obtain  all  the  details  of  the  angular  scattering  pattern.  To  illustrate  this 
point,  we  plotted  the  amplitudes  corresponding  to  x  =  4,  on  the  same  basis  as  in  Fig.  5. 
Figure  6a  corresponds  to  m  =  1.315  and  may  be  compared  with  a  similar  case  shown  by  van 
dc  Hulst.*  The  characteristic  large  amplitudes  at  small  scattering  angles,  where  =  Sj,  and 
the  rapid  decrease  in  amplitude  with  increasing  9  are  clearly  evident.  For  this  relative  size, 
a  computation  with  ^9  — 10°  is  necessary  for  a  graphical  interpolation  of  the  amplitude 
with  respect  to  #.  Figure  6b,  for  m  =  1.315  —  0.1370/,  shows  that  the  effect  of  absorption  in 
this  case  is  to  decrease  the  magnitude  of  the  amplitude  considerably  at  all  angles,  contrary  to 
the  effect  shown  on  Fig.  5.  Note  also  the  reversal  of  the  looping  in  52(^)  near  180°,  as  com¬ 
pared  with  Fig.  6a,  At  180°,  of  course,  we  always  have  5,(180°)  =  — Si(180°).  Increasing  the 
absorption  still  further  (Fig.  6r)  results  in  still  smaller  amplitudes,  which,  near  the  forward 
direction,  have  much  smaller  imaginary  parts.  Note  again  the  peculiar  looping  near  the  back¬ 
ward  directions. 

A  general  feature  appearing  in  Fig.  5  and  Figs.  to  6c  is  that  the  radius  vector  from 
the  origin  to  the  forward  scatterpoint  seems  to  rotate  clockwise  with  increasing  absorption 
index. 

The  amplitudes  for  larger  particles  (x  =  6)  are  partly  illustrated  in  Figs,  la  and  lb  (pp. 
22-23) .  These  smooth  curves  have  been  fitted  to  points  computed  with  A0  =  2°.5.  Larger  in¬ 
tervals  would  not  be  sufficient  to  yield  the  correct  behavior  of  the  function  by  smooth  curve 
fitting,  especially  in  the  backward  hemisphere.  The  number  of  loops  has  increased,  and  in  the 
absorbing  case.  Fig.  lb,  there  seems  to  be  a  cusp  in  5,(^)  near  9  =  90°.  Van  de  Hulst’s  remark 
that  one  full  turn  in  such  curves  takes  an  interval  in  0  of  the  order  360°/x  is  roughly  cor¬ 
roborated  in  Fig.  la  but  not  in  Fig.  lb.  Needless  to  say,  the  amplitude  curves  become  in- 


^Sce  p.  MS  of  Ref.  1. 
tSte  p.  23S  of  Ref.  1. 
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(O) 

Fig.  6 — Representation  in  the  complex  plane  of  the  amplitudes  Sx  and  5,  for  spheres  with  x  =  4. 
The  curves  for  the  same  indices  as  in  Fig.  ^  are  shown  separately  in  (a),  {b),  and  (r),  covering  the 
full  range  from  0°  to  180^  at  intervals  of  in  B. 


creasingly  complicated  for  larger  values  of  x,  and  the  intervals  in  9  must  be  made  even 
smaller  to  get  an  adequate  description. 

Examples  of  numerical  values  of  the  real  and  imaginary  parts  of  5^  and  5,  for  a  given 
m  and  0  =  0^  and  180°,  as  a  function  of  x,  arc  reproduced  in  the  tables  in  Appendix  C.  These 
show  the  variability  of  the  amplitude  at  these  angles  with  changes  in  the  size  of  the  sphere 
(or  changes  in  the  frequency,  if  the  sphere  size  is  constant),  for  various  indices  of  refraction. 

In  general,  the  amplitude  curves  derived  from  the  exact  theory  show  that  it  should  be 
relatively  easy  to  replace  the  Mie  series  by  analytical  approximations  even  in  the  absorbing 
case,  provided  0  is  kept  small.  Otherwise,  one  should  not  expect  too  much  from  simple  in- 
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Fig.  S—confintttd 
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Fic.  7 — ^The  amplitudes  Si  and  ^3  for  spheres  with  x  =  6,  computed  at  intervab  of  in 
(4)  Index  of  refraction  nt  =  1.315,  range  30°  ^  180°.  The  value  5(0°)  =  35.566  +  13669i 

is  not  shown  as  it  falls  outside  the  frame,  (b)  Index  m  =  1.315  —  0.4298/,  range  22°. 5  ^  ^  ^  180°. 
The  value  5(0°)  =  21.459  —  1.4633/  is  not  shown. 


terior  field  approximations  in  the  intermediate  frequenqr  range.  The  scattered  radiation  field 
depends  critically  on  the  interaction  of  reflected,  diffracted,  transmitted,  and  surface  waves. 
One  can  hardly  expect  to  reproduce  these  interactions  by  ad  hoc  separations  of  the  field  with¬ 
out  reaching  a  self-defeating  degree  of  complexity  in  the  mathematical  expressions. 

VARIATION  OF  THE  INTENSITY  AND  POLARIZATION  OF  THE  SCATTERED  LIGHT 

WITH  SCATTERING  ANGLE 

The  normalized  intensities  defined  by  (7)  through  (10)  are  the  most  important  param¬ 
eters  in  problems  of  radiative  transfer  for  polydispersed  media,  provided  independent  or  in- 
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Fig.  7 — continued 


coherent  scattering  can  be  assumed.  When  the  sphere  is  illuminated  by  unpolarized  radiation, 
the  scattered  radiation  shows  partial  linear  polarization  with  the  maximum  electric  field  oscih 
lation  either  perpendicular  or  parallel  to  the  scattering  plane.  For  convenience,  we  may  define 
the  polarization  P  in  this  case  by  the  expression 

/,((?) -/,(e) 

so  that  positive  and  negative  values  correspond  respectively  to  the  two  cases  mentioned  above. 
The  denominator  in  (14)  always  defines  the  total  intensity  of  the  scattered  radiation,  regard- 
less  of  its  state  of  polarization. 

The  remaining  parameters  /]|(^)  and  /4(^)  are  related  to  those  elements  of  the  scattering 
matrix  that  determine  the  Stokes  parameters  defining  the  plane  of  polarization  and  degree  of 
ellipticity  of  the  scattered  light  when  the  incident  light  is  polarized.^  These  parameters  are 
important  in  problems  with  multiple  scattering  even  when  the  incident  light  is  neutral.  They 
will  not  be  discussed  here,  since  one  has  to  specify  various  types  of  poiarizati<m  in  the  inci¬ 
dent  light  in  order  to  understand  their  effect  in  changing  the  polarization. 

Figures  8a  to  8r  arc  logarithmic  plots  of  /\(9)  and  as  a  function  of  scattering 
angle  S,  all  for  x  =  3  but  progressively  increasing  absorption  index.  In  all  three  cases  it  is 
seen  that  when  the  diameter  of  the  sphere  is  approximately  equal  to  one  wavelength,  the 
scattering  pattern  is  highly  asymmetrical,  and  the  polarization  deviates  considerably  from 


^See  pp.  40 f,  of  Ref.  1. 
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{O) 

Fic.  8 — The  logarithm  of  the  intensity  functions  /,(0)  (dashed  line)  and  i^(9)  (full 
line)  for  spheres  with  x  =  5,  plotted  at  intervals  of  in  6^  and  joined  by  straight-line 
segments.  The  curves  for  the  three  indices  m  =  1.315,  1.315  —  0.1370/,  and  1.315  — 
0.4298/  are  shown  separately  in  (4),  (b),  and  (r),  respectively. 


the  Rayleigh  case.  Light  scattered  at  angles  0^  <  0  <[  40°  exceeds  the  back-scattered  light  at 
angles  l40°  <  0  <  180°  by  two  to  four  orders  of  magnitude.  It  is  interesting  to  note  that 
this  ratio  does  not  depend  on  the  absorption  properties  in  a  simple  manner;  for  example, 
the  back-scattered  intensity  for  f'  =  0.4298  (Fig.  8c)  is  less  than  that  for  r'  =  0  (Fig.  84), 
but  much  larger  than  that  for  /  — 0.1370  (Fig.  8^).  A  broad  secondary  intensity  maximum 
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{C) 

Fig.  8 — continued 


between  90^  and  120^  is  evident  in  all  cases  shown,  but  it  is  less  pronounced  in  the  most 
absorbing  case. 

More  interesting  are  the  effects  of  absorption  on  the  polarization.  The  sign  of  polariza¬ 
tion  for  0  <  90°  is  the  same  for  all  three  cases,  but  in  the  backward  hemisphere  it  depends 
on  the  magnitude  of  the  absorption  index.  Furthermore,  the  degree  of  polarization  changes 
considerably  with  absorption.  In  the  regions  of  maximum  intensity,  especially  in  the  forward 
hemisphere,  the  polarization  seems  to  increase  directly  with  absorption,  as  shown  in  Table  2. 


RESULTS 
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Table  2 

P  hr  x  =  3,  f=:  1.315 


;  / 

0 

1 - 

0 

n - 

0.0143 

1  0.137 

1  .... 

1  0.4298 

40* 

0.098 

0.105 

0.159 

0.281 

60* 

0.041 

0.053 

.  -  1 

0.166 

0.460 

100® 

0,563 

0.593 

O.8O3 

0.995 

The  situation  for  larger  spheres  corresponding  to  x=  10  is  shown  in  Figs.  9a  and  9b. 
The  case  with  a  pure  dielectric,  m=  1.29,  in  Fig.  9a,  shows  a  number  of  secondary  maxima 
as  already  indicated  by  the  detailed  work  of  Penndorf  and  Goldberg^“*  and  reproduced  in 
Ref.  7.  Note  that  the  polarization  is  very  weak  and  negative  in  the  diffraction  maximum  at 
0°  <  <1  <  15®,  and  that  it  changes  sign  several  times  at  other  angles.  When  a  large  absorp* 


{O) 

Fig.  9 — ^The  logarithm  of  the  intensities  (dashed  line)  and  i,  (full  line)  for  spheres  with  x  =  10, 
plotted  at  intervals  of  2^.5  in  $.  The  curves  for  m  =  1.29  and  1.29  --  0.4720i  are  shown  separately 
in  (4)  and  (b),  respectively. 
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(b) 

Fig.  9 — continued 


tion  is  introduced  as  in  Fig.  9^,  the  forward  intensity  peak  remains  the  same;  however,  its 
polarization  becomes  positive.  Furthermore,  the  polarization  in  this  case  remains  positive  at 
all  angles  up  to  160®.  The  scattered  intensity  field  for  120°  <  <  180°  is  very  flat  com¬ 

pared  with  the  dielectric  case,  because  of  the  suppression  of  the  components  arising  from 
multiple  internal  reflections. 

As  mentioned  before,  with  the  addition  of  absorption,  polarization  in  the  forward 
hemisphere  increases  in  the  positive  direction.  This  effect  becomes  more  pronounced  with 
increasing  size  and  is  most  noticeable  in  the  regions  of  maximum  intensity.  G>mpare  the 
polarization  values  in  Table  3  with  those  given  for  x  =  3,  with  reference  \o  the  scattering 
diagrams  in  Figs.  Sa,  Sb,  8c,  9a,  and  9b.  In  the  region  40°  to  60°  and  for  x=  3,  the  polar¬ 
ization  in  the  most  absorbing  case  is  3  to  10  times  that  of  a  dielectric,  while  in  the  same 
region  but  for  x=  10  the  increase  in  polarization  is  by  a  factor  of  18  to  23. 


RESULTS 
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Table  3 

P  for  x=  10,  f  =  1,29 


/ 

$ 

0 

0.0472 

0.4720 

2y 

0.023 

0,031 

0.122 

4y 

0.018 

0.089 

0.457 

62^5 

0.042 

0.170 

0.762 

85® 

0.137 

0.382 

0.964 

0 

o 

0.409 

0.625 

0.701 

A  physical  explanation  of  this  in  terms  of  ray  optics  may  be  postulated  as  follows:  As 
mentioned  above,  the  intensity  and  polarization  of  the  scattered  light  is  the  sum  of  several 
streams  of  light  emerging  in  a  given  direction.  The  Fresnel  reflection  coefficients  for  the 
intensities  of  the  two  polarized  components  /j  and  /„  respectively,  after  one  reflection,  may 
be  written  in  the  form 


R  I.  (?  -  0*  + 

R  (qr  -  1)»  +  (fO* 


(15) 

(16) 


where  ^  =  sin  r/sin  /,  and  t  and  r  are  the  inclination  angles  of  the  incident  and  refracted 
rays,  respectively.  (Both  of  these  expressions  reduce  to  (12)  when  q=V,  that  is,  at  normal 
incidence.)  Keeping  q  constant,  the  examples  in  Table  4  demonstrate  the  changes  in  |  R,  |’, 
I R,  I’,  and  P  with  absorption,  obtained  from  (15)  and  (16).  These  values  show  that  in  gen¬ 
eral  external  reflection  of  the  unpolarized  incident  stream  should  add  positive  polarization, 
so  to  speak,  to  the  scattered  field.  This  effect  should  be  more  pronounced  at  the  smaller 
angles  of  scattering  where  the  magnitude  of  the  reflection  coefficient  is  large.  The  example 


Table  4 


P  from  external  reflection  for  r=  t.29 


f  = 

:0.5(#  40*) 

f  = 

=  0.8(«  =  80*) 

/ 

■qb 

■QB 

P 

■EB 

■HIM 

P 

0 

0.192 

0.047 

0.61 

0.055 

t).00025 

0.99 

0.472 

0.247 

0.066 

0.58 

0.101 

0.054 

0.50 

% 
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further  shows  that  at  ^  40^,  the  externally  reflected  rays,  in  both  the  nonabsorbing  and  the 

absorbing  cases,  should  be  positively  polarized  in  the  same  amount  (0.6l  and  0.58,  respec¬ 
tively).  However,  the  actual  polarizations  of  the  scattered  light  at  the  same  angle  (see  Figs. 
9a  and  91^)  are  — 0.07  for  the  nonabsorbing  and  0.53  for  the  absorbing  case.  The  negative 
polarization  in  the  nonabsorbing  case  may  be  thought  of  as  the  resultant  of  externally 
reflected  rays  plus  those  transmitted  and  reflected  internally  into  the  same  direction.  When 
the  latter  are  suppressed  by  absorption,  the  positive  polarization  of  the  externally  reflected 
ray  dominates  the  scattered  field.  In  the  present  example,  /  =  0.472  corresponds  to  an  ex¬ 
ponential  absorption  coefficient  y,  given  by 

47r/  5.89 

For  example,  in  visible  light,  "rays'*  penetrating  the  sphere  are  practically  extinguished  along 
a  path  of  a  micron  or  so  in  the  interior  \of  the  medium  because  of  the  large  value  of  the 
coefficient  y.  This  line  of  reasoning  also  explains  why  the  effect  of  absorption  on  the 
polarization  is  not  so  noticeable  in  the  case  of  smaller  spheres,  where  interior  paths  are  short 
and  geometric  optics  does  not  apply. 

Finally,  Fig.  10  illustrates  some  of  these  points  more  clearly  in  the  case  of  metallic 
spheres.  The  total  intensity,  +  /j,  and  the  polarization  P  are  shown  separately  for  an  iron 
sphere  illuminated  by  X  5890  A  radiation  (m  =  1.51  —  1.63/)  and  a  dielectric  sphere  with 
m~  1.525,  both  having  a  diameter  of  about  2.2  wavelengths.  The  six  secondary  maxima  in 
intensity  observed  in  the  dielectric  case  are  completely  suppressed  in  the  case  of  an  iron  sphere. 
The  polarization  curves  in  the  inset  of  Fig.  10  show  several  maxima  and  minima  in  both 
cases;  however,  the  iron  sphere  scatters  positively  polarized  light  in  all  directions,  whereas 
the  dielectric  sphere  displays  mostly  negative  polarizations.  These  differences,  particularly  in 
the  polarization,  are  extremely  important,  for  example,  in  the  investigation  of  the  nature 
of  interplanetary  particles  responsible  for  the  zodiacal  light. 
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*  =  7.0  ..4 
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-m  =1.525 


7  MYc™ 


20  40  60  80  100  (  20  140  160  180^ 

e 


•m  =1.525 


I""* 

'-/n « 


m  «  1.51  -  1.65^* 


100  120  140  160  180^ 

ScQtttring  onglt,  9 

Fig.  10 — The  logarithm  of  the  total  intensity  for  spheres  with  x  =  7,  plotted  at  intervals 

of  5®  in  $.  The  curves  for  an  iron  sphere  (m  =  1.51  —  1.63/)  and  a  dielectric  sphere  (i»  =  1.525) 
are  shown  on  the  same  diagram  for  comparison.  The  inset  shows  the  degree  and  sign  of  partial 
linear  polarization  of  the  light  scattered  by  each  sj^iere  in  various  directions  when  illuminated  by 
unpolarized  light. 


SUMMARY  AND  COMMENTS 


All  the  functions  related  to  the  scattering  of  electromagnetic  waves  on  homogeneous 
spheres  with  a  complex  index  of  refraction  can  be  easily  computed  to  a  high  degree  of  accu¬ 
racy  by  modern  machines.  Thus,  the  effects  of  absorption  on  the  mechanism  of  scattering, 
especially  in  the  intermediate  frequency  range,  can  be  exactly  assessed  without  the  use  of 
simplifying  assumptions  based  on  physical  and  geometrical  optics.  Such  approximations  have 
proved  successful  in  the  high-  and  low-frequency  cases,  but  they  generally  fail  in  the  middle 
range  except  under  severe  limitations  on  the  physical  properties  of  the  substance  in  the 
sphere.  It  is  hoped  that  the  present  results,  which  are  based  on  the  exact  theory,  besides 
revealing  hitherto  unknown  details  of  the  scattering  process,  will  provide  reliable  criteria 
for  checking  the  range  of  validity  of  any  extensions  of  approximate  theories  into  the  inter¬ 
mediate  range. 

But  there  is  a  still  more  important  application  of  this  type  of  study.  The  intermediate 
range,  especially  when  the  physical  properties  in  the  sphere  differ  considerably  from  those 
of  the  surrounding  medium,  is  precisely  the  range  of  interest  in  a  class  of  problems  involving 
the  absorption  and  the  oiffuse  transmission  and  reflection  of  sunlight  by  planetary  atmos¬ 
pheres,  including  our  own.  In  fact,  as  is  known  from  the  terrestrial  atmosphere,  the  maximum 
linear  size  of  solid  or  liquid  particles  that  can  remain  suspended  for  appreciable  lengths  of 
time  is  limited  to  a  few  tens  of  wavelengths  of  visible  light.  It  is  further  known  that  a  small 
fraction  by  mass  of  such  **large’'  particles  mixed  in  a  gas  is  sufficient  to  alter  radically  the 
character  of  the  transmitted  and  reflected  light.  Similar  remarks  apply  to  interplanetary  and 
interstellar  dust  particles.  Of  course,  the  assumption  of  spherical  particles  may  not  be  correct, 
but  it  is  useful  to  understand  completely  their  scattering  properties  before  investigating  more 
complicated  shapes. 

In  any  case,  a  complete  quantitative-theoretical  understanding  of  the  spectrum,  intensity, 
and  polarization  of  the  light  diffusely  reflected  by  sunlit  planetary  atnx>spberes  is  essential 
at  this  time  in  order  to  interpret  correctly  the  observational  data  obtainable  by  means  of 
instrumented  space  probes  sent  to  the  vicinity  of  the  nearer  planets.  In  'the  examples  pre- 
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seated,  we  have  shown  that,  in  addition  to  changes  in  the  size  of  the  spheres,  small  changes 
in  the  absorption  index  have  large  effects  on  the  intensity  and  polarization  of  the  scattered 
light,  particularly  at  large  angles.  In  the  planetary  problem,  a  consideration  of  models 
assuming  a  number  of  particles  that  are  identical  in  size  and  nature  is  hardly  sufficient. 
Such  a  situation  cannot  exist  either  in  planetary  atmospheres  or  in  interplanetary  and  inter- 
stellar  space.  The  ultimate  purpose  of  these  calculations  is  to  derive  the  scattering  properties 
of  a  small  sample  volume  of  sp^ce  containiilg  a  representative  distribution  in  the  size  and 
the  index  of  refraction  of  the  particles.  The  resultant  values  of  the  elements  of  the  scattering 
matrix  pertaining  to  the  aggregate  will  be  summations  of  those  pertaining  to  the  individual 
particles.  This  summation  cannot  be  accomplished  without  a  knowledge  of  the  scattering 
properties  of  individual  particles.  Both  theory  and  existing  measurements  show  that  the  scat¬ 
tering  properties  of  a  polydispersed  medium  tend  to  be  smooth  functions  of  the  scattering 
angle.  The  latter,  once  determined,  can  be  treated  analytically  to  separate  the  angular  de¬ 
pendence  from  the  size  dependence  in  order  to  solve  the  corresponding  equation  of  radiative 
transfer. 
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The  recursion  formulas  used  for  the  angle-dependent  functions  are  as  follows: 

^  “  „  lri  < 

r„(e)  =-.  cos  e  [^„(6f)  -  ]  -  (2«  -1)  sin*  e  7r„.,(0)  + 

with 

^o{^)=0,  To(fl)=0, 

7r,(e)  =  l,  T,(0)  =  COS0, 

^2(0)  ^  3  COS  0,  T.^(0)  =  3(cos*  9  —  sin’  9) . 

The  scattering  functions  were  then  computed  by  the  summation  of  the  absolutely  con- 
verging  series  indicated  by  the  expressions  (1)  to  (4)  in  the  text. 


APPENDIX  B 


RANGE  OF  COMPUTATIONS 


The  various  indices  of  refraction  and  size  ranges  covered  by  our  computations  are  summar* 
ized  in  Table  5,  where  r  represents  the  real  part  and  /  the  imaginary  part  of  the  refractive 
index,  written  w  =  f  —  //,  and  and  x,  represent  the  lower  and  upper  limits  of  the  rela¬ 
tive  size,  respectively. 


Table  5 


Refractive  Index 

Limits  of  Size 

— 

Wavelength 

(m) 

r 

Substance 

I.IU 

0.1831 

0.5 

10.0 

water 

11.50 

1.212 

0.0601 

0.5 

10.0 

water 

10.00 

0 

0.5 

water 

2.25 

1.29 

0.0472 

0.5 

water 

8.13 

0.0645 

0.5 

0.4720 

5.25 

0 

0.5 

— 

water 

1.61 

1.315 

0.0143 

0.5 

water 

5.30 

0.1370 

0.5 

water 

6.05 

0.4298 

0.5 

mSM 

water 

15.00 

1.34 

0 

0.3 

70.0 

water 

0.43 

1.353 

0.0059 

0.3 

1 

20.0 

water 

3.90 

1.44 

0 

0.3 

■■ 

0.4000 

0.5 

■■ 

water 

16.6 

0 

0.3 

■Si 

1.525 

0.0050 

0.5 

0.0100 

0.3 

0.0682 

0.3 

23.0 

water 

.  3.07 

1.28 

1.37 

0.1 

7.0 

iron  ^ 

0.441 

1.51 

1.63 

0.1 

7.0 

iron 

0.389 

1.70 

1.84 

0.1 

i _ 

7.0 

iron 

0.668 
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The  complete  angular  range  from  0°  to  180°  was  covered  in  steps  of  5°  or  2^^°,  except 
in  the  forward  area  between  0°  and  15°,  where  the  steps  were  made  smaller  to  get  the  details 
of  the  aureole. 


APPENDIX  C 


SELECTED  COMPLEX  AMPLITUDE  VALUES 


A  TABULATION  of  the  valucs  of  the  complex  amplitude  at  0  =  0°  and  180°  is  presented  in 
Table  6  for  selected  indices  of  refraction,  as  a  function  of  x  in  steps  of  0.)  in  x.  The  six 
significant  figures  of  the  IBM  printouts  have  been  rounded  off  to  four  in  this  table. 
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